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An  oven  has  been  designed  and  built  for  the  production  of  metal  micro- 
crystals via  nucleation  and  growth  of  the  metal  vapor  in  an  inert  gas 
atmosphere.  It  has  a flowing  argon,  double  orifice  sampling  system  producing 
a supersonic  free  jet  mixture  which  crosses  a 40  kev  electron  beam.  Debye - 
Scherrer  diffraction  patterns  are  obtained  from  the  metal  cluster  samples 
which  range  in  size  from  40  to  95  % in  diameter  with  estimated  concentrations 
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of  10  "-10  cm  J at  the  electron  beam  location.  The  average  cluster  size 
produced  in  the  oven  for  all  three  metals  studied  correlated  well  with  the 
product  of  oven  pressure  times  metal  evaporation  temperature  p0Tom-  As 
Po^om  increases,  average  size  increases  and  cluster  concentration  decreases, 
in  qualitative  agreement  with  other  metal  evaporation  research  and  with 
nucleatlon  studies  of  vapor-inert  gas  expansions  in  supersonic  nozzles. 
Analysis  of  the  diffraction  patterns  reveal  changes  in  crystal  structure 
from  that  of  the  bulk  in  the  neighborhood  of  50-60  X diameter  (2,000  to 
4,000  atoms  per  cluster).  The  differences  are  most  pronounced  in  indium 
which  changes  from  tetragonal  to  face  centered  cubic  as  the  size  decreases. 
Cluster  temperatures  at  the  electron  beam  location  have  been  estimated  from 
Debye-Waller  factors,  from  extrapolation  to  bulk  structure,  and  from  calcu- 
lated values  of  the  argon  gas  temperature  in  the  free  jet.  The  three  methods 
yield  consistent  results  and  the  hotter  the  clusters  are  prior  to  the  free 
jet  expansion,  the  greater  is  the  temperature  difference  between  the  clusters 
and  Che  argon  gas  when  the  mixture  becomes  collisionless  (i.e.  no  further 
change  in  the  temperature  of  gas  or  clusters). 
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INTRODUCTION 


>¥ 

The  objective  of  the  work  presented  here  is  to  determine  the  differences 
in  structure  of  small  metal  microcrystals  from  that  of  the  bulk  phase  as  the 
average  size  is  diminished.  The  size  at  which  changes  occur,  as  well  as  the 
nature  of  these  changes,  is  of  importance.  The  metal  is  vaporized  in 
an  oven  in  the  presence  of  argon,  an  inert  or  ^carrier*  gas  which  initiates 
the  metal  nucleation  and  growth  in  the  gas  phase,  near  the  evaporating 
surface.  The  newly  formed  microcrystals  or  clusters  are  entrained 
into  a subsonic  free  jet  and  transported  through  a double  orifice  sampling 
arrangement.  The  portion  of  the  argon-metal  cluster  mixture  that  passes 
through  the  second  orifice  or  nozzle  expands  as  a supersonic  free  jet  and  is 
crossed  by  a 40  kev  electron  beam  1 to  2 nozzle  diameters  downstream 
where  diffraction  patterns  are  taken  on  glass  photographic  plates. 

The  metal  cluster  source  is  designed  to  initiate  the  nuc 
"stagnation"  or  pre-expansion  chamber  instead  of  the  usual  technique  of 
expanding  an  unsaturated  vapor  through  an  orifice  or  nozzle.  In  the  latter 
case  the  phase  change  occurs  in  the  supersonic  expansion.  This  choice  was 
greatly  influenced  by  the  reported  difficulty  in  getting  metals  to  nucleate 
in  adiabatic  expansions.  For  example,  mercury  was  reported  as  not  clustering 

appreciably  in  nozzle  beam-type  expansions*1  and  in  large  laval  nozzles  with  cooling 

2 

rates  hundreds  to  thousands  of  times  slower  than  those  of  the  small  free  jets 
(i.e.  10'*  °C/sec  in  contrast  to  10^  to  10^  °C/sec  in  nozzle  beams).  One  might 
predict  such  behavior,  at  least  qualitatively,  on  the  basis  of  the  classical 

nucleation  theory  where  the  rate  of  formation  of  clusters  large  enough  to 

3 

initiate  the  condensed  phase  is  given  by 


J - K exp(-AG*/kT)  . (1) 


2 


Here  J Is  the  number  of  "critical"  size  clusters  formed  cm  sec  , k is 

Boltzmann's  constant,  T is  temperature  and  K is  a large  pre-exponential  factor 

dependent  on  the  thermodynamic  properties  of  the  supersaturated  state.  The  energy 

* 3 2 2 

barrier  for  this  Arrhenius- type  rate  equation  is  AG  /T  in  S,  where  a is 

surface  tension  and  the  saturation  ratio,  S = (PV/P®)T>  is  the  ratio  of  the 
vapor  pressure  to  the  equilibrium  vapor  pressure  at  the  same  temperature. 

Because  metals  have  surface  tensions  tens  to  hundreds  of  times  higher  than 
most  gases,  it  requires  enormous  supersaturations  to  lower  the  energy  barrier 

down  to  values  encountered  with  gases.  This  can  be  accomplished  by  vaporization 
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into  an  inert,  cool  gas  atmosphere  ’ . This  technique  also  greatly  reduces  the 
requirements  of  high  oven  power  and  large  vacuum  pumping  capacity. 

OVEN  DESIGN 

The  oven  cluster  source  of  Fig.  1 represents  a design  evolution  over 
several  previous  metal  sources  used  in  our  laboratory  in  which  the  gas  flow 
has  been  minimized.  In  addition,  the  necessary  heater  power  is  now  low 
enough  so  that  a highly  regulated  power  supply  can  be  used,  requiring  only 
minor  adjustment  of  the  electron  beam  with  no  jitter  due  to  AC  fields.  This 
design  also  optimizes  the  cluster  concentration  at  the  crossed-beam  inter- 
section or  scattering  volume,  providing  several  orders  of  magnitude  increase 
over  previous  oven  conf lgurations . 

The  source  is  built  to  provide  continuous  x - y motion  during  operation, 
relative  to  the  electron  beam.  The  x-motion  for  the  oven  as  a unit  is  pro- 
vided by  the  screws  Sj  of  Fig.  1.  Between  experiments  the  oven  can  be  dis- 
assembled and  the  distance  between  orifices  and  N2  changed  by  adjusting  the 
screws  S2  shown  also  In  Fig.  1.  During  an  experiment  the  heater  or  metal 
source,  argon  supply  Jet,  and  the  thermocouple  can  be  continuously  moved  as  a unit 
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relative  to  the  orifices,  or  nozzles  N^  and  N^.  Varying  this  distance 
changes  the  flow  time  from  the  heater  to  the  first  nozzle  Nj.  See  Table  I 
for  dimensions  and  characteristic  operating  conditions. 

The  oven  is  operated  with  either  a tungsten  filament,  or  a boron  nitride 

crucible  with  a coiled  filament  heater.  Either  configuration  can  be  viewed 
using  a window-prism  arrangement.  This  is  used  to  inspect  the  operation 
visually  as  well  as  providing  access  for  an  optical  pyrometer.  The 
argon  supply  and  the  pumping  rate  for  the  annular  region  between  the  orifices 
can  be  controlled  so  as  to  vary  the  pressure  in  the  oven  p^  or,  holding  pQ 
constant,  to  alter  the  jet  velocity  which  controls  the  gas  mixing  time.  The 
oven  temperature  is  varied  by  adjustment  of  the  filament  power  supply  and 
measured  with  the  optical  pyrometer.  The  voltage  drop  across  the  filament 
and  the  amperage  are  recorded  so  that  the  power  dissipated  in  the  oven  is 
known.  The  thermocouple  which  senses  the  gas  temperature  before  it  enters 
the  first  nozzle  has  a radiation  shield  over  it  to  avoid  direct  radiant 
heating  of  the  Junction  by  the  hot  filament. 

NUCLEATION  AND  GROWTH  IN  THE  OVEN  SOURCE 

Most  of  the  experiments  reported  here  were  conducted  using  a triple  strand, 
coiled  filament  across  the  copper  heater  leads  with  the  metal  forming  a 
spherical  droplet  of  about  0.3  cm  diameter  near  the  center  of  the  filament. 

A few  experiments  with  lead  were  performed  in  a 0.5  cm  diameter 
x 1 cm  deep  crucible.  The  nucleatlon  and  growth  zone  for  the  crucible  is 
one-dimensional,  while  for  the  filament  the  zone  is  one-dimens ional  only  if 
the  zone  location  and  thickness  is  small  compared  to  0.3  cm.  Otherwise,  it 
will  be  three-dimensional.  If,  in  either  configuration  the  flow  in 
the  nucleation  zone  is  turbulent  it  would  be  most  difficult  to  conceive  a 
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meaningful  theoretical  model  for  the  process.  Turbulence  alone  without  the 
complications  of  heat  transfer,  mass  transport,  and  nucleation  is  difficult 
enough  to  describe  theoretically.  Estimates  of  the  Reynolds  number  Re,  for 
the  argon  flowing  out  of  supply  line  near  the  heater,  range  from  10  to  50, 
which  is  well  within  the  laminar  flow  regime  (Recall  Re  = pvf/p,  where  p is 
gas  density,  v is  velocity,  l is  a characteristic  length  and  p.  is  the  dynamic 
viscosity).  Since  the  argon  velocity  over  the  filament  or  the  crucible  source 
is  expected  to  be  significantly  lower  than  that  at  the  argon  supply,  plus 
increased  viscosity  at  the  higher  temperatures,  the  Reynolds  number  near  the 
evaporating  metal  surface  is  even  lower  than  the  above  values.  There  has  been 

research  to  describe  theoretically  the  nucleation  in  a laminar  boundary 
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layer  ; and  in  no  flow  situations,  i.e.  stagnant  films  or  layers  ’ ; and 

9 10 

diffusion  cloud  chambers  ’ . The  case  of  the  stagnant  film  nucleation  has 

g 

many  features  in  common  with  adiabatic  nozzle  expansions  . If  the  oven 

11  12 

produces  convection  whether  natural  or  forced  , the  flow  geometry  can 
become  complicated,  again  making  theoretical  analysis  difficult. 

Due  to  the  primary  objective  of  this  work,  i.e.  structure  of  small  micro- 
crystals formed  in  the  gas  phase,  in  addition  to  the  above  mentioned  diffi- 
culties, a detailed  description  and  analysis  of  the  metal  nucleation  and  growth 

process  has  not  been  undertaken.  (Metal  nucleation  and  growth  for  bismuth, 

13 

lead,  and  iron  has  been  studied  in  shock  tube  experiments  .)  Certain 
gross  features  however  are  worth  pointing  out. 

1)  The  Increase  in  average  size  with  oven  temperature  (or  metal  vapor  density) 

5 14 

reported  by  others  ’ has  also  been  observed  in  this  work.  There  is  one 
Interesting  difference  In  the  size  determination  however.  Electron 
diffraction  Is  used  in  the  gas  phase  which  probes  the  microcrystal  size  that 
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results  from  nucleatlon  and  growth  only.  Electron  micrographs  may  or  may  not 

measure  the  primary  microcrystal  size  in  samples  if  coalescence  of  the  micro- 

2 

crystals  has  occurred.  Also,  our  argon  pressure  is  a factor  of  10  to  10  lower 
thus  reducing  the  tendency  for  the  microclusters  to  coalesce^. 

2)  The  increase  of  cluster  size  with  ambient  argon  pressure  follows  the  trend 
observed  in  other  evaporation  studies^.  This  feature  can  be  qualitatively 
explained  since  an  increase  in  argon  pressure  inhibits  the  metal  vapor  mass 
transport  shifting  the  nucleation  and  growth  regions  toward  the  higher  vapor 
density  region  thus  enhancing  growth.  This  is  similar  to  the  explanation 
pr  to  account  for  increased  particle  size  with  inert  gas  molecular 

.e.  increased  collision  cross-section'*  * . 

le  results  of  items  1)  and  2)  above  have  been  summarized  in  Fig.  2a. 

The  variation  of  microcrystal  "average"  size  with  both  the  metal  vapor  temp- 
erature T (and  therefore  its  pressure  and  density)  and  argon  pressure  pQ  is 
plotted  for  the  three  metals  investigated.  This  correlation  is  intended  only 
as  a qualitative  indication  since  the  actual  process  is  governed  by  nucleation 
and  growth  requiring  greater  detailed  knowledge  of  many  more  thermodynamic  and 
flow  properties  than  just  pQ  and  T^. 

4)  An  estimate  of  cluster  production  rate  and  concentration  at  the  electron 
beam  location  has  been  made.  A few  representative  experiments  for  each  metal 
are  used,  again  only  to  indicate  the  approximate  magnitude  and  trend  of  the 
cluster  production.  An  upper  limit  estimate  of  the  cluster  flux  or  production 
rate  is  derived  assuming  that  all  the  heater  power  is  taken  up  by  the  heat 
of  vaporization  for  the  evaporating  metal,  neglecting  losses  due  to 
radiation  and  to  heat  conduction  to  the  argon.  Dividing  the  resultant  mass 
flux  by  the  "average"  cluster  size  obtained  from  the  diffraction  pattern 
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ring  widths  gives  the  cluster  production  rate  Nc^  in  clusters-sec  , which 
is  shown  in  Fig.  2b. 

Using  the  known  spreading  of  the  metal  clusters  in  the  supersonic  free 

jet  flow  from  N,,  (2.5°  half  angle),  the  location  of  the  electron  beam 

(x  ,/D  „ ~ 1.5).  the  terminal  Mach  number  which  occurs  at  or  prior  to 
eb  nZ 

x/D^  = 1.5  (Hj,  = 2.5),  and  the  cluster  conservation  relation  = NcjAv  with 

A the  jet  cross-sectional  area  and  v its  velocity,  the  concentration  at 

the  scattering  volume  defined  by  the  intersection  of  the  free  jet  and  the 

electron  beam  can  be  determined.  If  the  cluster  flux  used  is  the  maximum  value 

then  the  cluster  concentration  will  be  the  maximum  or  upper  limit  also.  These 

calculated  values  of  Nc ^ have  been  reduced  by  an  order  of  magnitude  as  an 

estimate  to  account  for:  a)  only  a fraction  of  the  electrical  power  is  actually 

used  for  evaporation,  b)  some  of  the  clusters  will  be  transported  to  the  walls 

of  the  oven  source,  and  c)  a fraction  of  the  total  argon-metal  cluster  mixture 

flow  is  pumped  out  of  the  annulus  between  and  Nj.  Thus  an  estimated  cluster 

concentration  in  the  scattering  or  diffraction  volume  is  obtained  and  is 

suown  in  Fig.  2b  with  the  right-hand  ordinate. 

As  with  Fig.  2a  these  results  are  only  qualitative.  However,  one  obtains 

12  13  -3 

an  estimate  of  concentration  in  the  range  of  10  to  10  clusters  cm  and 
finds  a decrease  in  Nc ^ with  the  product  of  oven  pressure  and  oven  temperature, 
PoTotn»for  all  three  metals.  Thus  a general  picture  of  the  nucleation  and 
growth  of  metal  clusters  in  the  oven  source  emerges.  The  argon-metal  vapor 
mixture  cools  rapidly  with  distance  from  the  evaporating  liquid  metal  surface. 

In  spite  of  the  decreasing  vapor  density  with  distance  from  the  evaporating 
surface  the  metal  vapor  supersaturates  increasingly  with  distance  until  it 
reaches  a value  high  enough  to  cause  a nucleation  and  growth  process  that  is 


7 


rapid  enough  to  deplete  the  vapor  density  and  thus  the  supersaturation,  shutting 
off  the  nucleation  process.  The  remaining  metal  vapor  condenses  onto  the 
clusters  just  formed.  From  the  point  where  the  nucleation  rate  approaches  its 
maximum, the  clusters  formed  may  begin  agglomeration  or  coalescence  as  well  as 
growth  by  vapor  deposition.  From  all  indications,  theoretical  as  well  as 
experimental,  the  nucleation  and  grow  zone  or  boundary  layer  is  only  a few 
tenths  of  a centimeter  thick.  As  the  metal  surface  temperature  increases  so 
does  its  local  vapor  pressure  and  density.  The  supersaturation  profile  moves 
in  toward  the  surface  causing  the  maximum  nucleation  zone  to  move  with  it. 

The  higher  density  of  the  metal  vapor  here  causes  the  clusters  just  nucleated 
to  grow  faster  thereby  shutting  down  the  nucleation  sooner.  The  net  result 
then  is  a lower  concentration  of  clusters  and  since  there  is  more  vapor 
present  they  grow  to  larger  sizes.  Thus,  increasing  the  metal  vapor  pressure 
leads  to  the  formation  of  fewer  clusters  of  larger  size. 

Increasing  the  inert  gas  pressure  (or  inert  gas  collision  cross-section, 
i.e.  heavier  molecular  weight)  enhances  the  heat  transfer  while  inhibiting 
the  mass  transport.  This  also  shifts  the  supersaturation  curve  in  toward  the 
evaporating  surface  with  the  same  results  as  above,  a smaller  number  of  larger 
clusters  formed.  The  structure  of  the  nucleation  and  growth  in  a one- 
dimensional, stagnant  film,  thermal  boundary  layer  has  been  treated  by 
7 8 

several  authors  ’ . The  similarities  to  condensation  in  supersonic,  adiabatic 

g 

nozzle  flows  are  considerable  . The  trends  noted  above  are  identical  to  those 

observed  in  moist  air  expansion  in  nozzles  using  pressure  and  light  scattering 

measurements:  i.e.  increasing  the  condensable  species  partial  pressure  results 

17  18 

in  fewer  clusters  formed  which  then  grow  to  larger  sizes  ' 
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SAMPLING  THE  ARGON-METAL  CLUSTER  MIXTURE:  GASDYNAM1CS 

The  argon  supply  jet  issues  from  a copper  tube  of  0.15  cm  inside  diameter 

which  has  been  closed  down  at  the  end  to  0.1  cm  making  a small  converging 

section.  The  pressure  at  the  inlet  to  the  15  cm  long  supply  tube  is 

typically  1.6  torr  when  p^  = 1 torr.  Thus  the  pressure  at  the  tube  exit, 

though  unknown,  is  close  to  the  1 torr  value  in  thf  oven.  The  jet  leaving 

the  tube  is  thus  subsonic  and  mixes  with  the  ambient  argon  around  the  jet, 

setting  into  motion  a convective  flow  in  the  remainder  of  the  oven.  In 

particular  the  metal  vapor  which  nucleates  and  grows  in  a layer  near  the 

filament  heater  (or  within  the  crucible  as  the  case  may  be)  is  transported 

by  this  jet-induced  circulation  to  the  edge  of  the  jet  where  it  is  entrained 

as  the  jet  mixes  with  the  convection  flow.  The  mixture  is  pumrad  out  through 

nozzle  N^  to  a pressure  of  p^  = 0.6  torr  so  the  flow  through  N^  is  also 

subsonic.  If  it  were  isentropic  then  p^/p^  = 0.6/1  gives  a Mach  number 

M = 0.9.  Beyond  the  orifice  N^  the  flow  spreads  out  and  some  of  it  is 

pumped  off  with  a mechanical  vacuum  pump  of  140  liter  min  ^ capacity.  The 

mass  flow  to  this  pump  can  be  varied  with  a valve.  Then,  in  conjunction 

with  a valve  controlling  the  supply  of  argon  at  the  inlet,  the  oven  pressure 

p^  and  total  argon  mass  flow  rate  m^  can  be  controlled  independently.  This 

provides  control  of  p^,  an  essential  parameter  for  metal  nucleation,  as 

well  as  m which  determines  the  sampling  or  flow  time  from  the  heater  to 
a 

the  nozzle  N^.  The  intent  then  is  to  draw  in  the  mixture  using  the  mechanical 
pump  and  sample  it  with  nozzle  N2-  Under  the  actual  choice  of  operating 
conditions  used,  however,  the  greater  fraction  of  the  total  mass  flow  goes 
through  N^,  i.e.  typically  807.  is  seen  in  Table  I.  The  pressure  p^  is  not 
measured  directly  but  Inferred  from  the  background  pressure  p2  in  the  crossed 
beam  chamber  using  an  argon  flow  calibration  previously  made. 
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Using  known  mass  flow  rates  and  the  appropriate  dimens  ions  estimates 
of  the  gas  velocities  and  flow  times  in  the  various  chambers  have  been  made. 

Note  the  different  time  scales:  15  to  400  milliseconds  in  the  oven  chamber, 

10  to  100  microseconds  in  the  intermediate  pumping  chamber, and  2 microseconds 
from  to  the  electron  beam. 

The  dimensionless  ratio  of  viscous  forces  to  dynamic  forces  is  the  Reynolds 
number  Re  and  has  been  calculated  for  Table  I at  the  inlet  to  the  oven,  and 
at  nozzles  and  N2«  Low  values  of  Re  are  usually  associated  with  laminar 
flow  (e.g.  in  pipes  Re  < 2000).  However  even  though  the  flow  in  the  argon 
supply  tube  is  laminar  the  free  jet  may  be  mixing  turbulently  downstream  of 
the  filament  and  in  the  conical  region  of  the  intermediate  pumping  chamber. 

The  recirculation  region  in  the  oven  chamber  behind  the  argon  supply  jet 
as  mentioned  previously,  is  laminar  with  a Reynolds  number  much  less  than  20. 

The  ratio  of  the  argon  mean  free  path  to  a characteristic  dimension  is 
the  Knudsen  number, Kn  * \/l.  It  serves  to  delineate  the  regions  of  collision 
dominated  or  continuum  flow,  transition  flow,  and  free  molecular  or  collision- 
less flow  by  having  values  « 1 , sa  1 , and  » 1 respectively.  Thus  it  is 
seen  that  the  flow  in  the  oven  and  through  the  nozzle  Nj  is  continuum 
while  the  supersonic  flow  through  N2  is  approaching  the  transition 
regime.  There  are  two  important  features  which  derive  from  the  fact  that 
the  expansion  through  N2  is  approaching  a Knudsen  number  of  unity.  One  is 
that  the  x-position  in  the  free  jet  at  which  the  flow  becomes  collisionless 

is  very  close  to  the  nozzle  or  free  jet  orifice^.  The  position  at  which 

12 

collisions  effectively  terminate  is  x/Dn2  1 or  x • 0.075  cm  with  a typical 
value  of  terminal  Mach  number  M^,  - 2.5  and  a terminal  temperature  TT  » 117°K. 
Thus  placing  the  electron  beam  any  further  downstream  than  x/Dn2  “ 1 will 


10 

not  result  in  the  measurement  of  any  physical  change  in  the  clusters  due  to 
collisions,  such  as  colder  temperature  or  increased  size.  The  only  change 
with  increasing  x/D^  is  the  undesirable  trend  to  lower  cluster  concentration. 
The  second  feature  of  importance  is  that  due  to  the  relatively  small  number  of 
collisions  during  the  jet  expans  ion , the  jet  does  not  spread  out  as  much  as  a 
lower  Knudsen  number  flow.  The  metal  clusters  spread  through  only  a cone 
half  angle  of  2.5°  as  determined  by  collecting  a metal  film  deposit  at  two  dif- 
ferent locations  in  the  free  jet,  at  x = 0.15  cm  and  x « 50  cm.  This  is 
important  for  good  resolution  of  the  Debye- Scherrer  rings  in  the  diffraction 
patterns.  If  the  vertical  extent  of  the  free  jet  containing  the  metal 
clusters  becomes  an  appreciable  fraction  of  the  camera  length,  i.e.  the  dis- 
tance from  the  crossed  beams  to  the  film  plane,  the  peaks  would  be  additionally 
broadened  in  a manner  difficult  to  deconvolute. 

CLUSTER  COALESCENCE 

12  -3 

With  cluster  concentrations  in  the  range  of  10  cm  as  indicated  by 
Fig.  2b, and  characteristic  times  in  the  millisecond  range  or  longer,  coalescence 
is  predicted  theoretically  as  well  as  experimentally  observed  ’ . If 

coalescence  occurs  early  in  the  nucleation  and  growth  zone  and  the  clusters 
are  hot  enough  to  still  be  liquid  it  is  possible  that  they  form  larger 
homogeneous  liquid  droplets  that  later  solidify^.  The  diffraction  patterns 
from  these  microcrystals  should  be  characteristic  of  the  final  droplet  size. 

If  on  the  other  hand  the  clusters  cool,  solidify,  and  subsequently  agglomerate 
they  may  adhere  to  one  another  due  to  local  surface  attractive  forces  and 
exhibit  diffraction  patterns  characteristic  of  the  primary  microcrystals 
(i.e.  those  formed  via  nucleation  and  vapor  growth  only)  as  observed  in 
polycrystalline  thin  films.  Electron  micrographs  of  coalesced  particles  may 
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under  some  circumstances  reveal  the  primary  cluster  size  and  in  other  cases 
only  show  the  final  agglomerated  size. 

In  all  of  the  experiments  conducted  in  this  investigation,  whether  with 
lead,  bismuth,  or  indium,  in  which  cluster  size  is  determined,  there  is  a 
consistent  relationship  of  average  size  with  the  properties  which  most  directly 
affect  the  nucleation  and  growth  zone,  namely  metal  vapor  temperature 
and  argon  pressure.  On  the  other  hand,  a variation  by  a factor  of  two  in 
sampling  length  from  the  heater  to  nozzle  , causing  a factor  of  two  change 
in  the  longest  characteristic  flow  time,  produced  no  correlation  with  average 
cluster  size.  Although  these  correlations  were  not  pursued  in  an  extensive 
experimental  program, and  were  of  necessity  qualitative  due  to  the  demands  of 
our  primary  objectives,  it  is  nevertheless  tempting  to  conclude  that 
coalescence  did  not  significantly  affect  the  structural  features  of  the 
microclusters  which  are  presented  below. 

ANALYSIS  AND  RESULTS 

Diffraction  patterns  of  In,  Pb  and  Bi  clusters  are  recorded  on  photographic 
plates  (Kodak  Electron  Image  - 3^;  x 4"  plates)  under  the  conditions  described 
in  the  previous  sections.  Typical  microdensitometer  traces  are  illustrated 
in  Fig.  3 (In  and  Pb)  and  Fig.  4 (Bi).  Under  the  present  experimental  con- 
ditions, all  diffraction  patterns  show  relatively  sharp  Debye -Scher re r 
rings,  i.e.  no  obvious  liquid,  molecular  or  small  n-mer  (dimer,  trimer,  etc.) 
patterns  have  been  observed.  The  peak  positions  of  diffraction  patterns  can  be  con- 
sistently assigned  with  the  Miller  indices  (hk2)  of  the  corresponding  bulk 
crystals:  In  (D^h  - I4/mmm:  body  centered  tetragonal),  Pb  (0^  - Fm3tn:  face 

centered  cubic)  and  Bi  (D^  - R3m:  hexagonal).  In  the  case  of  In,  however, 
diffraction  patterns  change  their  contour  from  the  tetragonal  to  the  f.c.c. 
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as  the  average  size  of  cluster  becomes  small  as  illustrated  in  Fig.  3.  Once 
the  diffraction  peaks  are  indexed,  the  values  of  the  lattice  parameter  can 
be  calculated  from  the  following  equations. 


In  : 

x/LX  - 

{2(h2  + k2)/a2  + jf2/c2j* 

(2) 

Pb  : 

x/LX  = 

(h2  + k2  + /2)^/a 

(3) 

Bi  : 

x/LX  - 

f 2 2 2 2 23  ^ 

^(h  + k + hk)/3a  + L I c j 

(M 

2 2 

x - r - 3r  /8L 

(5) 

where  r is  the  radius  of  the  Debye -Scherrer  ring  indexed  as  (hki) , and  L\  is 
the  camera  constant  determined  by  a standard  thin  metal  film  sample  (L  = 53  cm). 
The  body  centered  tetragonal  indium  can  be  regarded  as  a distorted  f.c.c. 
structure  or  f.c.  tetragonal:  c/a  “■  1.075  in  a bulk  crystal  at  room  temperature. 
In  Eqn . (2)  for  In,  the  lattice  parameters  (a  and  c)  correspond  to  those  in 
the  distorted  f.c.c„  while  for  the  f.c.c.  indium  patterns  Eqn.  (3)  is  applied. 

For  Bi  known  as  a rhombohedral  structure,  a hexagonal  unit  cell  is  used  to 
determine  the  lattice  parameters. 

The  average  cluster  size,  D,  is  determined  using  the  relation: 

D - KX/0  cos  6 *“  L\/(W2  - W2)^  (0  « 1)  , (6) 

where  X is  the  electron  wavelength,  K is  the  shape  factor  (K  “■  1.0),  0 is  the 

Bragg  angle,  0 is  the  true  line  broadening  in  radians,  and  W and  W are  the 

o 

half-widths  of  the  lines  of  the  sample  and  the  reference  (Au  and  A l standard 
thin  films)  materials,  respectively.  For  the  estimation  of  D,  the  well  separated 
peaks  have  been  used:  (101)  for  In,  (111)  and  (220)  for  Pb,  and  (102)  for  Bi. 

The  lattice  parameters  and  average  cluster  diameter  D are  summarized  in 
Tables  2,3  and  4,  together  with  the  corresponding  experimental  conditions.  The 
lattice  parameters  and/or  their  ratio  are  plotted  In  Fig.  5 (In),  Fig.  6 (Pb), 
and  Fig.  7 (Bi)  as  functions  of  the  cluecer  site  (D)  and  the  number  of 
metal  atoms  per  cluster  (n). 
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Indium  Clusters 

A structure  change  mentioned  earlier  (Fig.  3)  in  In  microcrystals  can  be 
more  quantitatively  observed  in  Fig.  5.  The  variation  in  lattice  parameters  a 
and  c,  and  of  their  ratio,  c/a,  with  cluster  size,  shows  that  a bulk  f.c.  tetra- 
gonal structure  changes  into  f.c.c.  structure  near  D “ 50  X (n  3000).  This 
change  involves  a shrinkage  in  c and  an  expansion  in  the  basal  axis  a,  resulting 
in  the  unit  cell  volume  remaining  nearly  constant  (x  0.17.).  Another  interesting 
aspect  of  this  Figure  is  the  essentially  constant  values  in  a and  c,  and  thus 
also  in  the  ratio  c/a,  beyond  D > 65  X (n  ? 5000).  This  feature  permits  extra- 
polation of  the  lattice  parameter  data  to  large  sizes  and  application  of  known 

bulk  physical  properties  as  one  means  of  estimating  cluster  temperature.  The 

22 

thermal  expansion  coefficients  of  a bulk  In  crystal  are  known  for  a and  c 

® = 25.4  x 10'5  + 3.22  x 10’10  T2  and  a - 23.0  x 10-5  - 4.35  x 10'10  T2 

a c 

(80°K  ^ T < 300°K) . Using  the  extrapolated  values  of  a « 4.62  X and 
c * 4.89  X and  the  above  expansion  coefficients,  the  cluster  temperature  T^ 
is  estimated  to  be  264  - 46°K.  Similarly  T ^ can  be  obtained  from  the  extra- 
polation of  c/a  (1.058)  and  the  unit  cell  volume  (104.4  X~*) : T . = 362°K  and 

cl 

294°K,  respectively.  A simple  average  of  these  three  values  is  307°K. 

It  is  known  that  the  lattice  parameters  of  raicrocrystals  decrease  as  the 
23  24 

cluster  size  decreases  ’ . Such  behavior  for  rare  earth  metals  in  thin  film 

experiments  has  been  interpreted  in  terms  of  surface  tension  and  the  bulk 
24 

modulus  . It  should  be  noted,  however,  that  the  size  dependence  of  c/a  in 

the  present  result  can  not  be  explained  by  bulk  compression  due  to 

the  surface  tension,  since  in  bulk  In  crystals  c/a  increases  monotonically  with 

25 

an  increase  in  external  pressure  (e.g.  c/a  * 1.076  at  0 kbar  -•  c/a  - 1.096 
at  151  kbar),  while  the  present  results  show  that  c/a  decreases  as  the  cluster 
size  decreases  (i.e.,  as  the  pressure  due  to  the  surface  tension  Increases). 
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In  the  present  case  the  c/a  ratio  change  of  indium  may  arise  from  the  stability 
of  microcluster  structure.  For  small  clusters  surface  atoms  represent  a 
large  fraction  of  the  total  number  and  crystal  defects  such  as  vacancies  are 
more  likely  to  occur  in  the  surfaces  of  high  curvature.  These  defects  may 
play  an  important  role  with  regard  to  cluster  stability  dictating  a change 
from  tetragonality  to  cubic  indium. 

Lead  Clusters 

The  range  of  cluster  sizes  studied  for  Pb  goes  from  ~ 40  X in  diameter  to 

~ 80  X,  similar  to  the  case  for  In.  The  results  are  shown  in  Table  3 and  in 

Fig.  6.  In  contrast  to  In,  the  lattice  parameters  observed  in  this  range 

of  cluster  size  are  essentially  constant  within  an  experimental  error  of  «*  0.3%. 

The  figure  appears,  however,  to  indicate  a lattice  parameter  shrinkage 

below  D < 30  X (n  ^ 2000).  An  analysis  of  the  Bragg  peak  intensities  for 

the  smallest  clusters,  D “•  40-50  X,  shows  a vestige  of  an  amorphous 

solid;  the  detailed  result  will  be  published  elsewhere.  T . is  estimated  to 

c a 

be  212  K from  the  extrapolation  of  the  lattice  parameters,  as  has  been  made 

in  In  (see  Fig.  6),  using  a linear  thermal  expansion  coefficient  of 

a » 27.08  x 10  **(°C)  * (90  < T < 287°K)^  . From  a damping  of  Bragg  peak 

intensities,  i.e.  Debye-Waller  factor  , T ^ for  the  largest  cluster  size 

has  been  estimated  to  be  150°K:  T . - (212  + 150) /2  - 181°K. 

cZ 

For  the  lead  experiments  two  different  types  of  oven  sources  are  used 
to  evaporate  the  metal;  one  is  a tungsten  wire  coil  similar  to  those  used  for 
In  and  Bi,  and  the  other  is  a boron  nitride  crucible.  The  open  circles  in 
Fig.  6 correspond  to  data  taken  with  the  crucible  oven,  where  the  cluster 

sizes  are  relatively  small,  end  show  little  variation  in  size  with  a change 

in  p T . Although  no  systematic  study  of  nucleation  kinetics  has  been 
o ow 

carried  out,  it  certainly  indicates  that  the  size  distribution,  i.e.,  the 
nucleation  and  growth  processes,  is  sensitive  to  the  type  of  oven. 
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Bismuth  Clusters 

Average  cluster  sizes  in  the  range  D ^ 60-95  X have  been  studied.  7ne 
results  are  given  in  Table  4 and  Fig.  7.  The  lattice  parameters  for  this 
range  of  sizes  are  essentially  constant  within  an  estimated  error  of  0.37. 

(note  the  expanded  ordinate  in  Fig.  7 relative  to  those  in  Fig.  5 and  Fig.  6). 
The  apparent  decrease  in  c/a  ratio  with  decreasing  cluster  size  is  merely 
suggestive.  By  comparison  with  the  results  for  indium  and  lead,  such  a 
change,  if  any,  would  be  expected  to  occur  below  D ^ 50  X.  The  cluster  tem- 
perature T ^ to  be  166  - 20°K,  by  the  use  of  extrapolated  values  for  lattice 
parameter  a of  4.54  X and  for  c of  11.83  X and  bulk  thermal  expansion  coef- 
ficients27: cra  = 11.8  x 10‘6(°C)'1(193  < T < 373°K)  and  »c  = 17.7  x 10‘6(°C)'1 

(150  ^ T ^ 540°K) . When  we  use  the  c/a  ratio  of  2.606  for  the  estimation  of 

o 

Tc£,  it  gives  a relatively  low  value  of  80  K,  while  the  unit  cell  volume 
extrapolation  ( 211.17  X^)  provides  T ^ = 176°K.  A simple  average  of  the 
above  is  141°K. 

DISCUSSIONS  AND  CONCLUSIONS 

The  grain  size  dependence  of  the  lattice  parameters  for  In,  Pb  and  Bi 
microclusters  has  been  studied  in  the  range  of  40  ^ D ^ 95  X.  Above  D > 60  X, 
the  lattice  parameters  are  nearly  constant  and  approach  those  of  the  bulk 
crystalline  materials.  In  the  vicinity  of  D ““  50  X (2000  ~ 3000  atoms  per 
cluster),  the  lattice  parameters  begin  to  deviate  from  those  of  the  larger 
size  clusters.  A striking  example  is  the  case  of  In,  where  a transition 
from  tetragonal  to  cubic  has  been  observed. 
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Cluster  temperatures,  obtained  from  the  lattice  parameter  extrapolation 
in  the  previous  section,  are  307°K  (In),  181°K  (Pb)  and  141°K  (Bi),  respectively. 
It  is  interesting  to  compare  these  values  with  those  estimated  from  gasdynamic 
considerations,  where  it  is  assumed  that  the  free  jet  expansion  is  isentropic 
and  that  the  cluster  temperature  is  equal  to  the  carrier  gas  terminal  tempera- 
ture. Recall  that  the  free  jets  in  this  work  reach  their  terminal  state  at 
x^/D^  = 1,  that  is,  only  one  nozzle  diameter  downstream.  From  that  point  on 
the  flow  is  essentially  collisionless  and  the  gas  and  cluster  temperatures 
can  no  longer  change  except  possibly  via  radiation.  Upon  computing  typical 
mass  flow  rates  for  the  metals  and  for  the  argon,  and  using  their  respective 
solid  and  gaseous  heat  capacities,  it  is  found  that  the  argon  in  this  two 
phase  flow  has  more  than  50  times  the  heat  capacity  of  the  metal  microcrystals. 
Thus  computing  the  argon  terminal  temperature  TT, using  the  isentropic  equations, 
involves  a negligible  error  due  to  the  heat  exchange  between  the  cluster  and 
the  argon  carrier  gas. 

For  the  largest  clusters,  such  an  adiabatic  cooling  has  been  calculated 

and  is  listed  in  Table  5:  161-229°K  for  In,  117°K  for  Pb  and  146-150°K  for 

Bi,  respectively.  Comparing  these  to  the  cluster  temperatures  listed  above, 

it  is  possible  to  examine  the  temperature  lag  between  the  clusters  and  the 

expanding  argon  at  the  position  in  the  free  jet  where  the  electron  beam 

is  located,  x . /D  . » 1.5.  Recall  that  this  location  is  downstream  of  the 
eb  n2 

position  for  the  terminal  Mach  number  so  that  the  diffraction  experiment  is 
sampling  a collisionless  two  phase  flow.  Figure  8 presents  the  temperature 
difference  AT  - ^cA~^T  as  * ^unct*on  cluster  temperature  for  the 
expansions  that  contain  the  larger  microclusters  (i.e.  those  approaching 
the  bulk  structure).  It  is  seen  that  AT  Increases  as  the  temperature  of  the 
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clusters,  and  thus  also  the  cluster  stagnation  temperature  prior  to  expansion, 
increases.  * The  clusters  that  start  out  hotter  end  up  lagging  the  final  argon 
temperature  by  a greater  amount. 

For  the  clustering  of  pure  gas  expansions  in  free  jet  molecular  beams  it 
30 

has  been  reported  that  the  cluster  temperature  is  linearly  proportional  to 
the  well  depth  of  inter-atomic  (or  molecular)  potentials.  Note  in  Fig.  8 that 
the  cluster  temperature  increases  from  Bi  to  Pb  to  In.  So  also  do  their  heats 
of  sublimation.  Thus  the  observed  temperatures  of  metal  clusters  correlate 
almost  linearly  with  the  binding  energy  (i.e.  heat  of  sublimation).  Due  to 
the  nature  of  the  cluster  nucleation  and  growth  process  in  this  work  the 
correlation  of  cluster  temperature  involves  a somewhat  indirect  connection 
with  the  heat  of  sublimation.  Although  the  relationship  is  not  direct  the 
correlation  is  provocative. 

In  conclusion,  the  present  nozzle  assembly  for  generating  metal  micro- 
clusters is  suitable  for  electron  diffraction  studies  of  metal  cluster 
structures.  Its  compactness,  ease  of  handling  without  the  need  of  a cooling 
system,  and  attainment  of  high  temperatures  with  low  electric  power, 
have  a good  deal  of  practical  merit.  In  addition,  a predictable  method  for 
controlling  the  cluster  size  distribution  has  been  found.  This  holds  promise 
for  the  study  of  much  smaller  cluster  sizes. 
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Fig.  1. 


Fig.  2. 


Fig.  3. 


Fig.  4. 


FIGURE  CAPTIONS 

Schematic  diagram  of  a metal  cluster  generator  for  electron 
diffraction:  TC  - thermocouple,  W - prism  mirror  window,  H - heater 

for  metal  vaporization,  and  - sets  of  screws  for  configuration 
adjustments,  Ar  - argon  carrier  gas,  and  N 2 - flow  orifices 
or  nozzles;  see  text  for  details. 

a)  The  average  cluster  size  D,  and  b)  the  cluster  production  rate  ^ 

and  the  cluster  concentrat ion  at  the  crossed  beam-intersection  N . are 

CJt 

plotted  as  a function  of  the  product  of  argon  stagnation 

pressure  and  the  temperature  at  the  surface  of  the  evaporating 

metal.  The  general  trend  is  that  as  p T is  increased  clusters 

o om 

of  larger  diameter  and  lower  concentration  are  formed. 

The  microdensitometer  traces  of  electron  diffraction  photographs 
for  indium  and  lead  microclusters  show  tetragonal  In  of  an  average 
cluster  size  of  ~ 63  X in  diameter  and  is  compared  with  cubic 

(f.c.c)  In  of  D ^ 45  X.  A typical  f.c.c.  pattern  of  lead  is  shown 

with  an  average  size  of  ~ 45  X.  The  abscissa  S is  the  scattering 
parameter  S = (4Tf/X)sin  9 with  9 the  Bragg  angle. 

A typical  diffraction  pattern  of  bismuth  clusters  with  an  average 
cluster  size  of  ~ 60  X in  diameter  is  plotted  as  a function  of  S. 

The  relative  intensity  here  is  not  calibrated  from  the  photographic 
density.  A rhombohedral  Bi  structure  is  assigned  by  the  corre- 
sponding hexagonal  cell.  Some  of  the  Miller  indices  have  been 
omitted  here  due  to  space  limitation. 


^ 1 
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Fig.  5. 


Fig.  6. 


Fig.  7. 


Fig.  8. 


Lattice  parameters  and  their  ratio  for  the  indium  microclusters 
are  plotted  as  a function  of  the  number  of  atoms  per  cluster  (n) 
and  the  average  cluster  diameter  D.  Broken  lines  indicate  the 
corresponding  values  in  bulk  crystalline  materials  at  various 
temperatures.  Typical  experimental  error  bars  are  shown  by 
vertical  lines.  Deviation  from  bulk  structure  begins  as  the 
diameter  approaches  60  X. 

Lattice  parameters  of  lead  microclusters  are  shown  with  open 
circles  for  data  from  the  crucible  oven  and  closed  circles  for 
data  taken  with  the  coiled  tungsten  wire  oven. 

Lattice  parameter  ratios  of  bismuth  microclusters  are  presented 
as  a function  of  n and  D. 

Typical  cluster  temperatures  at  the  crossed-beam  location,  as 
determined  from  (Debye-Waller  factors)  and  from  unit  cell 
parameters  are  plotted  as  a function  of  the  difference  in 
temperature  between  the  clusters  and  the  argon  carrier  gas, 

“ - Tc/-TI- 
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TABLE  2 

INDIUM  CLUSTER  RESULTS 


In 

D(X)a) 

c(X)b) c) 

•*b) 

c/a 

p (torr)C ^ 
o 

T (°K)d) 
o 

T (°K)e) f) 
om 

Exposure 
time (sec)  ^ 

1 

81 

4.903 

4.642 

1.056 

0.91 

568 

1344 

2 

2 

81 

4.880 

4.606 

1.059 

0.91 

450 

1205 

10 

3 

68 

4.891 

4.627 

1.057 

0.91 

460 

1205 

7 

4 

63 

4.874 

4 . 644 

1.048 

0.91 

568 

1344 

5 

5 

63 

4.871 

4.617 

1.055 

0.50 

420 

1310 

10 

6 

62 

4.848 

4.608 

1.052 

0.75 

440 

1310 

10 

7 

52 

4.751 

4.663 

1.019 

0.67 

395 

1144 

30 

8 

51 

4.712 

4.712 

1.00 

0.67 

376 

1144 

25 

9 

50 

4.707 

4.707 

1.00 

0.71 

420 

1158 

20 

10 

44 

4.731 

4.731 

1.00 

0.63 

410 

1144 

45 

11 

44 

4.704 

4.704 

1.00 

0.71 

410 

1144 

40 

12 

42 

4.712 

4.712 

1.00 

0.71 

445 

1158 

30 

Bulkg) 

• 

4.945 

4.600 

1.075 

at  298°K 

a)  average  cluster  diameter;  estimated  error  range  of  5-15  X, 

b)  lattice  parameters  for  face-centered  tetragonal  and  face-centered  cubic 
structures;  the  experimental  errors  are  estimated  to  be  0.3-0,01,. 

c)  argon  carrier  gas  pressure  in  the  mixing  region;  see  text. 

d)  ambient  temperature  at  the  mixing  region,  measured  by  a Aiumel -Chrome  1 
thermocouple;  see  Fig.  1. 

e)  tungsten  wire  temperature  measured  by  an  optical  pyrometer. 

f)  exposure  time  of  electron  diffraction  photographs,  beam  current  ca.  4.6  pA. 

g)  bulk  crystalline  material  data,  taken  from  Ref.  22. 
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TABLE  3 


LEAD 

CLUSTER  RESULTS8'* 

Pb 

D(X) 

a(X) 

P 

(torr) 

o 

T (°K) 
o 

T (°K) 
om 

Exposure 

time(sec) 

1 

82 

4.939 

0.83 

352 

1266 

30 

2 

60 

4.938 

0.67 

423 

1333 

20 

3 

51 

4.939 

0.50 

410 

1205 

30 

4 

40 

4.935 

0.55 

376 

1205 

20 

* 

5 

51 

4.933 

1.0 

640 

1283 

20 

★ 

6 

51 

4.929 

0.98 

616 

1166 

20 

★ 

7 

49 

4.937 

0.98 

688 

1339 

23 

* 

8 

44 

4.937 

1.08 

688 

1286 

30 

* 

9 

40 

4.932 

1.0 

640 

1283 

15 

Bulkb) 

ao 

4.9505 

at  298°K 

a)  see  the 
spond  to 
for  the 

footnotes  of 
data  taken 
rest  of  the 

Table 
by  the 
plates 

2.  The  plate  numbers  with  * corre 
use  of  a boron  nitride  oven,  while 
a tungsten  wire  oven  was  used;  see 

b)  bulk  crystal,  taken  from  Ref.  28. 
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TABLE  4 

BISMUTH  CLUSTER  RESULTS*^ 


Bi 

D(X) 

a(&) 

c(X) 

c/a 

pQ(torr) 

T (CK) 
0 

T (°K) 
om 

Exposure 
time (sec) 

1 

94 

4.54 

11.83 

2.606 

0.83 

415 

- 

20 

2 

78 

4.54 

11.83 

2.606 

0.53 

352 

1366 

40 

3 

60 

4.55 

11.81 

2.594 

0.55 

340 

1177 

20 

Bulkb) 

00 

4.546 

11.862 

2.609 

at  298°K 

a)  see  the  footnotes  of  Table  2, 
those  for  a hexagonal  cell. 

The  lattice 

parameters  a and 

c are 

b)  bulk  crystal  data  taken  from  Ref.  29. 


TABLE  5 


ESTIMATION  OF  CLUSTER  TEMPERATURE  FROM  GASDYNAMICS 


Sample 

/ \a) 

Po(torr) 

To(°K)a) b) c) 

u c> 

*r 

TT(°K)d) 

ln-1 

0.91 

568 

2.11 

229 

In-2 

0.91 

450 

2.32 

161 

Pb-1 

0.83 

352 

2.46 

117 

Bi-l 

0.83 

415 

2.31 

150 

Bi-2 

0.53 

352 

2.06 

146 

a)  argon  pressure  inside  the  nozzle. 

b)  ambient  temperature  in  the  mixing  region;  see  Figure  1. 

c)  terminal  Mach  number;  cf  Table  1. 

d)  terminal  gas  temperature  which  is  a lower  limit  for  the 
cluster  temperature,  i.e.  for  a finite  heat  transfer 
rate  from  cluster  to  the  cooling  gas  the  cluster  temper- 
ature is  greater  than  the  terminal  gas  temperature, 
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